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A modified version of susceptibie-infected-recovered-susceptibie (SIRS) modei for tfie outbreaks 
of foot-and-moutli disease (FMD) is introduced. Tlie model is defined on small-world networks, 
and a ring vaccination programme is included. This model can be a theoretical explanation for 
the nonlocal interactions in epidemic spreading. Ring vaccination is capable of eradicating FMD 
provided that the probability of infection is high enough. Also an analytical approximation for this 
model is studied. 
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I. INTRODUCTION 

Currently, there are some outbreaks of foot-and-mouth disease (FMD) in some developing and European countries. 
FMD is very infectious for cattle and sheep that are essential source for human food. The outbreaks of FMD cause a 
high loss in the national income, because export markets are closed for a country once a disease appears in it. 

Generally animals do not have natural immunity against FMD. Also, FMD can be transmitted through contacts 
vifith infected animals, tools, food and aerosol. So once an infected animal is discovered, one can consider the whole 
farm is infected. There are some strategies to control an outbreak. Mass vaccination means vaccinating the whole 
population with a certain rate. But vaccination is expensive and a virus usually has different configurations that re- 
quire different types of vaccination. Also immunity due to vaccination often vanishes within some months. Therefore, 
one has to concentrate in vaccinating animals in contact to an infected case. This is called ring vaccination [1,2]. It 
follows the following procedures: 

(1) Once an infected case is discovered, the animals at the infected farm are slaughtered. 

(2) Animals are vaccinated within a ring with a certain radius around the infected farm. 

(3) A region around the infected farm is closed off. 

Ring vaccination always cost less than mass vaccination. Also it has been proven that ring vaccination is more 
suitable for FMD than mass vaccination [2]. 

Generally animals live in farms or herds forming clusters. Once a farm is infected, its nearest neighbours become 
at a high risk of infection. But this does not mean the interaction is just locally. In the 1987/1988 outbreak that 
occurred near Hannover, Germany [2], first two infected farms were discovered. A mass vaccination programme had 
been applied. After 3 months, another 2 infected farms were discovered at 18 Km away from the initial outbreak. 
Some days later, a third infected case was observed at a distance of 8 Km from the initial outbreak. This ensures the 
presence of nonlocal interaction in the spreading of an outbreak. 

Regular graphs can display the clustering property only. On the other hand random graphs [3] can show the 
nonlocal interactions without clustering. Small-world networks (SWN) [4,5] are shown to combine both local and 
nonlocal interactions. Also the concept of SWN was applied to many different systems. The results are often closer 
to the real systems than using regular lattices [6-8]. 

A coupled map lattice (CML) [9] is a dynamical system with discrete time and space steps and continuous states. 
It is used to model systems that consists of coupled elements. These systems often display spatioteniporal chaotic 
behaviour. The concept of CML is used to explain some aspects in many systems in different fields like biology, 
economics, etc. [9]. 

There are many mathematical models were suggested to study the epidemic spreading in a population [10-13]. 
According to the health of each individual, the population is classified to susceptible (S), infected (I) or recovered 
(R), In this work, we begin by studying a simple SI model on SWN. Then a model describes the epidemic spreading 
of FMD is suggested. It is an SIRS model [10] defined on SWN and a ring vaccination programme is considered. Also 
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an analytical approximated model for the epidemic spreading on SWN is generalized to the inhomogeneous case using 
some concepts of CML. 



II. SMALL-WORLD NETWORKS 



The concept of SWN [4,5] is proposed to describe social networks. The social networks arc characterized by two 
main properties: clustering and small-world effect. Small-world effect means the average distance between any two 
vertices is too short in compared with the size of the lattice. Only SWN is shown to satisfy the two properties 
together [5]. It is a connected 1-dimensional lattice of size L, with periodic boundary conditions. Some randomly 
chosen vertices are joined by some shortcuts to randomly selected other lattice sites. Here we consider shortcuts with 
length fc = 1. Let (j) be the average number of shortcuts per bond in the lattice. Hence for large L, the probability 
that two random vertices are connected by a shortcut is V' — 2k(f)/L. Naturally the critical concentration of this graph 
Pc is smaller than that for the ring p^ = l, and this was derived in [11]. 

Motivated by a simple susceptible-infected (SI) model, the spread of an epidemic in a population with random 
susceptibility is studied on SWN. This work generalizes that of Moore and Newman [11]. Consider the population 
occupy the vertices of a SWN of size L. Each individual is assumed to have a random susceptibility 1 — such 
that p{i) e (0, 1). The propagation rule is: 

(i) S-individual i having at least one infected neighbour (shortcuting neighbours are included, if exist) is infected 
if p(i) < q, where q is the infection rate. 

It is relevant to estimate the number of infected persons in this model as a fimction of time [14]. Setting q = 1 
and p{i) = 0V«=1,2,...,L]. Hence the number of infected individuals will grow initially as a sphere with surface 
r^t"^"^, where Fi = 2, r2 = 2%, Fs = 47r and so on. This is called the primary sphere. The probability of finding a 
shortcut is 2cf)rdt'^~^ per unit time. Once a shortcut is reached a secondary sphere is formed and so on. Hence the 
total number of infected cases is given by 



V{t) =tJ T 

Jo 



t 



1 + 2(f)V{t - r) 



dr. (1) 



Defining V = 2(j)V, t = t[2(l)Td{d — l)!]^/** and differentiating d times with respect to t one gets 



Its solution is 



= 1 + y. 



1=1 ^ ' 

In one dimension, one has Vif) = exp(f) — 1, in 2-dimensions V{t) = cosh(t) — 1. For t <1, the number of infected 
individuals grow as a power law f^/dl. While for t > 1, it grows exponentially. The transition occurs at i = 1 i.e. at 
t = 2(f)Td{d — 1)!. This result has an important effect in vaccination policies. It implies that vaccination should be 
administered as early as possible and with the highest possible ability to avoid reaching the exponential phase. Also 
immunizing individuals with shortcuts is more efiicient than immunizing ordinary individuals. 



III. NUMERICAL SIMULATIONS 



The susccptible-infected-recovered-susceptible (SIRS) model [10] is more closer to FMD. The transition between 
the states S, I and R occur according to the following rules: 

(i) (Infection): An S-individual having at least one infected neighbour becomes infected in the next time step with 
probability q. 

(ii) (Recovery): I- individuals are recovered by a rate (73. 
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(iii) (Losing immunity): The recovered individual change to the S-class by a rate g2- 



As a mean field approximation, the time rate of S, I and R are 




qSI. 



(3) 



(4) 




(5) 



These equations have two steady states: 



Si =N, h= 



Ri =0. 



(6) 




l2 = 



92 +93 



S2), R2 = — 12- 



(7) 



To get a positive value for L2; N must be greatcir than 5*2. Then the disease can be eradicated if < 1. 

This mean field approximation ignores the spatial structure of the lattice, and is valid for the case of global inter- 
acting model. But in reality, an outbreak spreads locally with some nonlocal interactions. To combine both local and 
nonlocal interactions, we define this model on SWN. A 1-dimensional lattice with periodic bomidary conditions is 
assumed. A fraction of 2(j)% of its individuals is assumed to interact with a third individual (randomly chosen through 
the lattice), in addition to the local interaction with the nearest neighbours. Also a ring vaccination programme is 
included to control the disease propagation. So the following procedure is added to that of the SIRS model: 

(Ring vaccination): With a probability cr, an infected case and its nearest neighbours is changed to the R-state. 

The shortcuting neighbours are not included in this step, because they are usually unknowns. 

The simulations wc;re carried out on a SWN with L = 5000 and (f> = 0.05. The shortcuting neighbours are fixed 
beforehand. The values of q2 = 0.2 and qs = 0.05 are fixed through the simulations. The value of q is varied q £ (0, 1) 
and the corresponding smallest a sufficient to eradicate the disease is calculated. The results are averaged over 10 
independent nms and given in Fig. (1). 

An analytical approximation for the model is given in Appendix A. 



Using SWN is better than using regular or random lattices in modelling the outbreaks of FMD. SWN describe both 
local and nonlocal interactions. 

For a simple SI model defined on SWN: Initially the number of infected individuals grows as a power law, then 
after a critical value of time, it grows exponentially. So vaccination should be adminstered as early as possible to 
avoid reaching the exponential growth. 

Ring vaccination is capable of eradicating FMD even for high infection probability. 



Here a model is presented that approximates the epidemic spread on SWN. This model has been used before in 
different contexts [12,13]. Now it is generalized to the inhomogeneous case. Consider n patches each one contains 
a certain number of individuals (say animals). In general these patches are not identical. Infection spreads from 
infected animals within the patch and due to those diffusing from other patches. Then the time rate of the number 
of infected individuals in the ith patch is given by: 



IV. CONCLUSIONS 



APPENDIX A: ANALYTICAL CALCULATIONS 
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(Al) 



The first term represents the infection within the patch with a rate Aj. The second represents the effect of other 
patches both nearby and far away at a rate /i^. The recovery rate is represented by 7^. Since the effect of other 
patches (second term in Eq. (Al)) is significantly smaller than the first one we expect that fii <C A^. Now the disease 
is eradicated if = Vi = 1, 2, ...,n. So the stability of this solution is studied. The system (Al) is a kind of CML, 
so we begin by presenting general stability results for CML. Then it is applied to Eq. (Al). 
Typically a 1-dimensional CML is given by 



= (1 - D)e 



D 



9* 

0+1 



or 



D 



eY' = {i-D)f{e]) + - 



where t = 1,2, ... and j = 1,2, n. Consider the following inhomogeneous steady state of the system 



= a,- 



Linearizing around this solution, the system (A3) becomes 

D " 



-f+i 



{1-D)f\aj)s' + 



/'(ai-i)4-i+/'K+i)4+i 



(A2) 



(A3) 



(A4) 



(A5) 



A useful result on the eigenvalues of the system (A5) is Gerschgorin theorem [15]. It states that the eigenvalues A of 
a square matrix [aij] satisfy |A — aul < X^^^j , hence 



|A| < 



(A6) 



Hence the steady state solution (A4) is stable if the following conditions are satisfied 



\{1-D)f{a,)\ 
For the homogeneous case: 



-^f (0:3+1) 



+ 



<lVi = l,2,...,n. 



(A7) 



(A8) 



Then the stability conditions reduce to |/'(a)| < 1. Condition (A7) shows how diffusion may stabilize the system, 
since if f'{oij) > 1 and f'{aj±i) < 1. Then Eq. (A7) may be still valid for enough large coupling. The stability 
conditions for the homogeneous solution (A8) for the system (A2) are 



27rfc ■ 

1 — cos( ) 

n 



1-D 

Generalizing to the 2-species CML, we get 



+ /'(") 



< lVfc = 0,l,...,n-l. 



r^r = (1 - D)fj + f (/j+r + fu) , vr = (1 - m + f (^j+i + au) ^ 



(A9) 



(AlO) 



where /j = f{xj,yj), = g(a;*, y*). Consider the homogeneous steady state x* — xo,yj = yo- Linearizing around it, 
by taking a;*- = a;o + e*-, yj = yo + Vj- Also define the doublet: 



x5 = 



(All) 



Then the linearized equations can be written in the following form 
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(A12) 



where 



Ak = {l- D) 



dl df 

dx dy 

dg_ dg^ 

dx dy 



at {xo,yo)- 



(A13) 



This system corresponds to the direct product matrix C = Aj I + Aj ® T, + Aj T," ^ , where 



1 ... 

1 



(A14) 



1 ... 

Thus the steady state of the system (AlO) is stable if all the eigenvalues Afc of the matrix A satisfy |Afe| < 1, where 

A = Aj+ Aj+ie-2'^''=/" + Aj_ie2'^''=/". (A15) 
The equations of a CML with delay [16] are given by 



D 



= (1 - D)f{x]) + - + , yf 1 = x*. 



(A16) 



Repeating the above analysis, we obtain that the steady states a;*- = yj = xo are stable if all the eigenvalues of the 
following 2x2 matrix B = Bj + B^+ie-^'^''^/" + Bj.ie^'^''^/" satisfy \Xj^k\ < 1, where 



Bi 



il-D)f'{x,) 
1 



, Bj±i 



#/'fe±i) 





The stability conditions for the inhomogeneous solution (A4) of the 2-species system (AlO) is: 



(A17) 



(1-D)( 



df 



dxj-i 



df 



dy, 

df 



f( 



df 



dxj + i 



df 



dvn 



<1, Vj = l,2,...,n. 



(A18) 



(A19) 



For the delayed system (A16), the stability conditions of (A4) is: 

{l-D)\f;\ + ^\f;^,\ + ^\f;_,\<l, Vj = l,2,...,n, 
where /j = df/dxj. 

From the above results, it is straightforward to prove the following proposition: 

Proposition 1: If local components are stable, then the homogeneous steady state solution of the correspond- 
ing CML is stable. 



Globally coupled CML are typically given by: 
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4+' = (i-^)/(^*) + ;^E/(4)- 

The inhomogeneous steady state is stable if 



D 



(1 - D) \f{aj)\ + ^ E < 1' = 1' 2' 



(A20) 



(A21) 



Continuous time CML are typically given by: 



(A22) 



The homogeneous steady state is given by 9j = Ooij — 1, 2, n and d6/dt = 0. Hence /(^o) + 5(^*0) + ff(^o) — 0, and 
it is stable if 



Re 



/'(^) + 25'(^)cos 



27rr 



<0, r = 0,1,2, ...,n- 1. 



Another form of continuous time CML is 



rif) 

= m) + [9{ej+,) + 9{e^-,)]h{e^). 

In this case the steady state is /(^o) + b(^o) + g{6o)]h{6Q) = 0, and it is stable if 



Re 



2/Tr 

f'{0) + 2h'{0)g{e) + 2h{e)g'{0) cos — 



< 0, r = 0, l,2,...,n- 1. 



(A23) 



(A24) 



(A25) 



Returning to Eq. (Al), the stability of the zero solution is determined by the eigenvalues of the following matrix 



A = 



at 111 iJ,i ... Hi 

M2 0:2 M2 • • • At2 

IJ.3 113 as ... /X3 

/^n • • • - 



, ai = Xi- 7i. 



(A26) 



For the homogeneous case where all the patches are identical, the parameters are independent of i. Hence the matrix 
A is circulant and the largest eigenvalue is A — 7 + /x. Therefore if A — 7 < and A — 7 + > 0, then the disease 
would have been eradicated locally, but the diffusion term may cause it to persist. This shows the importance of the 
diffusion term (the long range edges in the case of SWN) . 

Considering the more realistic inhomogeneous case where the patches are different. Also taking the natural assump- 
tion that /ij -C Aj, the matrix A is studied keeping only up to the second order terms in fj,. Thus the Routh-Hurwitz 
stability conditions are that all the following determinants are positive: 



Al = ai, A2 = 



ai 1 
a2 



ai 1 
03 02 ai 
as 



where 



at 1 ... 

as 02 ai . . . 

A _ as a4 as . . . 

Q ... an 



ai = - ^ aj, a2 = ^ aiUj - ^/ij, as = - ^ {akaj - HkfJ-j), 

a4 = ^ OliOll (Q!feQ!j - MfeMj), . . . , 



a„ = (-1)" 



The effect of ring vaccination is expected to reduce Aj not fii. 
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Figure Caption 

Fig. 1: The relation between tlie infection probability q and the smallest probability a sufficient to eradicate a disease. 
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